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ABSTRACT 

This  report  deals  with  the  prediction  of  non-equilihrium  infrared 
radiation  signatures  in  the  wakes  of  re-entry  vehicles.  The  first  section 
deals  with  the  prediction  of  non-L.TE  wake  molecular  radiation  for 
ablating  vehicles.  The  second  section  discusses  the  radiation  signatures 
which  can  result  from  the  injection  of  small  particles  into  the  wakes. 
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I.  INTRODUCTION 

The  general  subject  of  this  report  is  the  prediction  of  non- 
equilibrium infrared  radiation  signatures  in  the  wakes  of  re-entry  vehicles. 
The  subject  matter  presented  is  broken  into  two  distinct  parts.  The  first, 
presented  in  Section  II,  deals  specifically  with  the  prediction  of  non- 
equilibrium wake  radiative  histories  of  ablating  vehicles.  A vibrational 
and  chemical  kinetic  scheme  appropriate  for  vehicles  with  heat  shields 
composed  of  the  elements  carbon,  hydrogen  and  fluorine  is  presented. 
Detailed,  two-dimensional,  non- equilibrium  radiation  predictions  for 
both  laminar  and  turbulent  wakes  are  provided  for  a generic  re-entry  vehicle 
over  the  altitude  range  of  69  - 38  km.  The  relationship  between  the  pre- 
dicted radiative  behavior  and  the  non-equilibrium  vibrational  and  chemical 
kinetics  is  discussed. 

In  Section  III  the  radiative  behaviour  which  can  result  from  the 
injection  of  small  particles  into  the  wake  is  examined.  A model  is  devel- 
oped to  describe  the  interaction  between  the  particles  and  wake  molecules. 
The  model,  valid  in  the  limit  of  free  molecular  flow,  includes  the  pheno- 
mena of  gas/particle  collision  dynamics,  vaporization  and  radiative  losses. 
The  model  is  then  exercised  for  typical  particle  materials  and  sizes  in 
order  to  elucidate  the  importance  of  the  gas/particle  interaction  to  both 
spoofing  and  masking  applications. 
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U.  INFRARED  WAKE  RADIATION  SIGNATURES 
OF  ABLATING  RE-ENTRY  VEHICLES 


A.  Introduction 

In  recent  years  considerable  effort  has  been  directed  towards  the 

prediction  of  the  infrared  radiation  histories  of  the  wakes  of  ablating 

vehicles.  Such  predictions  were  initially  performed^  under  the  assumptions 

of  frozen  wake  chemistry  and  vibrational/translation  equilibration.  More 
2,  3 

recently  detailed  predictions  have  been  provided  which  include  the 
effects  of  non-equilibrium  vibrational  and  chemical  kinetics.  These  latter 
calculations  have  been  limited  to  laminar  wakes  and  did  not  include  the 
appropriate  vibrational  kinetics  for  HF,  a species  which  can  be  of  some 
importance  for  specific  heat  shield  combinations. 

This  section  is  concerned  with  the  details  of  the  techniques  used  in 
performing  such  calculations  and  the  generic  dependences  of  wake  infrared 
histories  on  such  parameters  as  altitude  and  ablation  products.  The  basic 
model  used  for  the  calculations  is  described  in  Sec.  B.  Specific  vibrational/ 
chemical/radiative  predictions  for  a "typical"  re-entry  vehicle  are  presented 
in  Sec.  C along  with  an  overview  of  kinetic  considerations. 

B.  Modeling 

A basic  computer  code  package  used  for  the  prediction  of  laminar 
wake  infrared  radiation  histories  has  been  developed  over  the  last  several 
years  and  has  recently  been  documented  in  some  detail.^  The  basic  input 
to  this  package  is  the  ablation  history  for  a specified  trajectory  of  the  re- 
entry vehicle  of  interest.  A series  of  computer  codes  are  then  used  to 
determine  the  flux  of  the  various  chemical  species  which  enter  the  wake, 
with  the  individual  species  concentrations  being  determined  through  use  of 
a partial  equilibrium  boundary  layer  model.  These  species  fluxes  are  then 
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used  as  input  for  a two-dimensional  axi- symmetric  wake  code  which 

incorporates  the  non- equilibrium  vibrational  and  chemical  kinetics.  This 

code  is  the  heart  of  the  computer  package  and  provides  as  output  the  space- 

time  vibrational /chemical  histories  of  all  wake  species.  This  output  may 

then  be  used  in  conjunction  with  a radiation  code  to  provide  wake  infrared 

radiation  histories  in  specified  bandpasses. 

3 

This  basic  computer  package  was  taken  as  the  starting  point  for 

the  present  calculations,  however  a few  modifications  were  made.  In 

terms  of  the  computer  programs  the  major  i odification  was  performed 

on  the  two-dimensional  wake  code.  Since  this  code  requires  the  integration 

of  a large  number  of  "stiff"  differential  equations  in  two  dimensional  space 

(i.  e.  the  chemi/vibro- relaxation  equations),  the  computer  calculations  are 

both  lengthy  and  costly.  The  predictor-corrector  integration  scheme  used 

4 5 

in  the  code  was  originally  derived  by  Lewis,  who  has  recently  developed  a 
somewhat  more  sophisticated  technique.  The  computer  code  has  now  been 
modified  to  include  this  new  integration  scheme  and  it  has  been  found  that 
the  program  running  time  has  been  reduced  by  a factor  of  1.  7 with  no  loss 
in  accuracy. 

The  second  major  modification  in  the  computer  package  involved  the 

vibrational/chemical  kinetic  scheme  used.  This  scheme  was  originally 

2 

developed  several  years  ago  and  included  all  important  chemical  reactions 
involving  the  species  H,  O,  N,  CO,  H^,  NO,  N^,  OH,  O^,  CO^  and  H^O. 

With  regard  to  vibrational  kinetics  NO  and  OH  were  taken  to  be  in  equilibrium 
in  order  to  minimize  the  number  of  reactions.  All  other  vibrational  reac- 
tions of  importance  were  included.  In  the  modeling  each  important  vibra- 
tional mode  of  a molecule  was  treated  as  an  independent  species  and  only 
those  vibrational  modes  which  had  significant  populations  at  temperatures 
of  2500°K,  the  approximate  peak  wake  temperature,  were  included.  As 
might  be  expected  this  involves  a large  number  of  vibrational  reactions, 
some  of  which  had  not  been  studied,  and  a number  of  the  rate  constants  (or 


the  temperature  dependence  of  rate  constants)  had  to  be  estimated. 
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This  kinetic  package  has  been  updated  under  the  present  program. 

Furthermore  the  vibrational  kinetics  of  HF  have  been  added  to  the  scheme 

inasmuch  as  this  species  can  be  the  dominant  wake  radiator  in  the  3-4 |im 

wavelength  region.  The  update  was  an  indepth  one,  utilizing  the  recent 

6 7 

reviews  by  Lewis  and  Trainor  and  by  Taylor  as  well  as  a large  number 
of  relevam  measurements  published  in  the  last  few  years.  This  overall 
review  resulted  in  some  fifty  changes,  additions  or  deletions  in  the  kinetic 
package.  The  updated  kinetic  scheme  used  may  be  found  in  Appendix  A. 
There  are  of  course  still  a number  of  rate  constants  which  required 
estimation,  however  it  is  felt  that  for  the  most  part  the  dominant  reaction 
mechanisms  are  well  defined. 

There  is  one  point  which  should  be  made  relative  to  the  chemical 
reactions.  Since  each  vibrational  mode  is  treated  as  a distinct  species  the 
individual  chemical  reactions  of  excited  states  must  also  be  specified.  This 
turns  out  to  be  an  advantage  since  it  allows  a convenient  method  for  incor- 
poration of  chemi- excitation  effects  such  as  enhancement  of  reaction  rate 
constants  due  to  vibrational  excitation  of  the  reactants  and  the  creation  of 
reaction  products  in  non-eqviilibrium  state  distributions.  The  only  chemi- 

excitation  reactions  of  importance  in  the  present  kinetic  scheme  are  the 

♦ 

binary  processes  for  production  of  HF 

H + F^  -►  HF*  + F (1) 

F + -»  HF*  + H (2) 

Both  of  these  reactions  can  produce  HF  in  highly  excited  vibrational 
levels,  up  to  V = 6,  and  a very  large  kinetic  scheme  would  be  required  in 
order  to  include  all  reactions  of  importance  between  these  vibrational 

g 

levels.  In  the  present  work  it  has  been  assumed  that  reactions  (1)  and  (2) 
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produce  HF  only  in  the  first  vibrational  state.  Although  this  is  a poor 
approximation,  it  turns  out  that  for  the  cases  considered  chemical 
production  of  HF  in  the  wake  is  not  the  dominant  source  of  HF,  i.  e. , it 
originates  in  the  equilibrium  boundary  layer.  Care  should  be  taken  how- 
ever when  applying  this  kinetic  scheme  to  other  systems. 

C.  Predictions 

1 . Initial  conditions 

Wake  calculations  have  been  performed  over  an  altitude  range  of 
69  - 38  km  for  a generic  medium-sized  slender  vehicle.  Re-entry  velocity 
was  taken  to  be  approximately  7 km/sec  with  initial  wake  velocities  relative 
to  the  body  of  approximately  6 km/sec.  All  wake  calculations  were  initiated 
at  a temperature  of  2500°K  with  initial  temperature,  velocity  and  species 
profiles  specified  as  radially  constant  across  the  wake.  Past  experience 
has  shown  that  such  top-hat  profiles  rapidly  decay  to  the  standard 
Ga ussian -like  profiles  in  the  near  wake.  Wake  edge  conditions  were  as 

Q 

specified  by  the  standard  atmosphere  tables^. 

In  order  to  include  a large  variety  of  chemical  species  in  the  wake 
the  vehicle  was  taken  to  be  composed  of  a teflon  nose  with  a carbon 
phenolic  frustrum.  The  initial  near  wake  mole  fractions  (moles/original 
mole)  used  in  the  calculations  are  shown  in  Table  I.  These  are  deduced 
by  integrating  over  the  predicted  partial  equilibrium  boundary  layer  species 
profiles. 

For  the  case  considered,  fluorine -containing  species  represent  the 
dominant  contaminants  at  higher  altitudes  since  the  teflon  nose  is  heated 
to  ablation  temperatures  more  rapidly  than  the  phenolic  afterbody.  Thus 
at  the  higher  altitudes  most  of  the  hydrogen  is  tied  up  in  the  species  HF  and 
the  concentration  of  water  vapor  is  minimal.  As  the  altitude  decreases 
the  afterbody  heats,  the  level  of  carbon  phenolic  ablation  increases,  and 
relatively  large  concentrations  of  H2O,  and  OH  are  found  in  the  wake. 

It  should  be  noted  that  there  are  significant  concentrations  of  free  radicals 
such  as  O,  H,  F,  OH,  etc.  at  all  altitudes  and  thus  there  is  a strong  poten- 


tial for  chemical  reaction  in  the  wake. 
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SPECIES 

H = 69  km 

H = 6 1 km 

H = 53 

km 

H = 46 

km 

H = 38  km 

F 

1.10  X lO"^ 

1.03  X 

10-^ 

3.78  X 

10-’ 

6 . 76  X 

10-^ 

3.95  X lO""^ 

H 

1.04x  10"^ 

3.00  X 

10-^ 

9.82  X 

10-’ 

3.89  X 

10-^ 

7.08  X lO’^ 

O 

8.07  X 10“^ 

7.61  X 

io‘^ 

6.45  X 

lO"^ 

4.01  X 

10-^ 

2.24  j(  lO"^ 

CO 

6.52  X lO"^ 

1.27  X 

10-^ 

1.96  X 

10-' 

5.27  X 

10-^ 

9.68  X 10"^ 

^2 

1.00  X lO"^ 

1 .00  X 

10-  = 

1.00  X 

10-' 

1 .00  X 

10-' 

1.00  X 10“^ 

HF 

4.59  X lO"^ 

1.29  X 

10-^ 

2.08  X 

10-^ 

2.58  X 

10-^ 

2.42  X 10  ^ 

«2 

1 .17  X lO"^ 

1 . 1 6 X 

10-^ 

2.33  X 

Ifl-' 

6.87  X 

10-' 

2.78  X lO'^ 

^2 

7.71  X 10"^ 

7.60  X 

10-' 

7.48  X 

10-' 

7.00  X 

10-' 

6.41  X lO"^ 

NO 

7.88  X 10"^ 

8.97  X 

10-^ 

9.67  X 

10-' 

9.12  X 

10-' 

7.40  X lO"^ 

OH 

1.36  X lO"'* 

4.70  X 

lo-* 

2.74  X 

10-' 

1 .17  X 

10-^ 

1.38  X lO"^ 

°2 

1.58  X lO"^ 

1.53  X 

10“^ 

1 .49  X 

10-' 

1 .20  X 

10-' 

9.35  X 10"^ 

CO2 

2.95  X lO"^ 

5.42  X 

10'^ 

8.23  X 

10-' 

1 . 87  X 

10-^ 

2.16x10 

H,0 

6.23  X lO"^ 

7.20  X 

10-' 

1.23  X 

10-' 

1 .43  X 

10-^ 

2.53  X 10"^ 

j 
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Non- equilibrium  chemical/vibration  wake  predictions  have  been 
made  at  the  five  altitudes  listed  in  Table  1.  The  calculations  at  the  highest 
four  altitudes  were  performed  for  the  case  of  a laminar  wake  with  the 
standard  molecular  viscosities  applied  to  the  diffusion  equations.  The 
calculation  at  46  km  was  carried  to  a shorter  distance  downstream  than 
those  at  higher  altitudes  in  order  to  account  for  the  fact  that  the  length  of 
laminar  wake,  i.e.  distance  behind  the  vehicle  prior  to  turbulent  transition, 
decreases  with  decreasing  altitude. 

At  the  lower  altitude  considered,  38  km,  turbulent  transition  was 
assumed  to  occur  directly  behind  the  body  and  the  two-dimensional  wake 
code  was  modified  in  order  to  allow  the  prediction  of  turbulent  wake  charac- 
teristics. This  was  accomplished  by  providing  an  option  to  use  a turbulent 
rather  than  laminar  viscosity  in  the  calculation.  The  turbulent  eddy  vis- 
cosity applied  was  that  developed  by  Donaldson  and  Gray,^^  defined  by 
the  relationship 


€ 


1 

0.  08  r 


(3) 


where  the  subscripts  e and  c correspond  to  wake  edge  and  centerline 
properties  respectively.  In  the  above  equation  p is  density,  U w’ake  velocity, 
r radius  and 


«l/2 


u + u 

c e 


Lastly  r is  that  value  of  wake  radius  where 


U (r')  = (U  - U ) e"^  + U 
c e < 


(4) 


(5) 
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2.  Chemical  kinetics 

Examples  of  the  effect  of  chemistry  on  wake  species  concentrations 
are  shown  in  Figs.  1-5.  Wake  centerline  histories  of  the  mole  fraction  of 
H^O  at  the  five  altituies  considered  are  displayed  in  Fig.  1.  At  the  highest 
altitudes  shown  wake  chemical  kinetic  times  are  typically  slow  because  of 
the  low  wake  densities,  thus,  significant  concentration  changes  due  to 
chemical  reaction  are  not  to  be  expected.  There  is,  however,  a noticeable 
increase  in  the  H^O  concentration  at  the  altitudes  of  61  and  69  km,  mostly 
resulting  from  the  reaction 


OH  + OH-^H^O  + O (6) 

This  increase  occurs  because  H^O  is  a minor  ablation  product  species  at 
these  altitudes  and  therefore  a small  production  rate  can  produce  a significant 
concentration  change.  Indeed  it  can  be  seen  that  at  the  lower  altitudes  of 
53  and  46  km  where  the  initial  H2O  concentration  is  higher  there  is  little 
change  in  the  H^O  concentration.  This  is  partially  due  to  the  fact  that  the 
centerline  wake  temperature  is  ~2300  - 2500°K  over  this  length  of  wake. 

Off  wake  axis,  where  the  temperature  is  lower,  there  is  some  production  of 
H2O  at  these  altitudes. 

The  H^O  mole  fraction  prediction  at  38  km  is  quite  different  than 
those  for  higher  altitudes.  The  drop  in  H^O  concentration  is  a result  of  rapid 
dilution  resulting  from  the  fast  growth  of  the  turbulent  wake.  Although  it  j 

cannot  be  deduced  from  the  centerline  history  shown,  there  has  been  a sig- 
nificant conversion  of  H and  OH  to  H2O  at  this  altitude.  In  other  words  the 
total  flux  of  H^O  in  the  wake  has  increased;  although  the  local  mole  fraction 
is  smaller,  the  total  wake  area  has  increased  in  compensation.  ^ 

Centerline  histories  of  the  mole  fractions  of  CO  and  CO2  for 
altitudes  of  61  and  46  km  are  shown  in  Fig.  2.  At  the  highest  altitude  the 
[CO]  and  [ CO2]  mole  fractions  both  decrease  slowly  due  to  laminar  diffusion. 

Since  the  total  wake  flux  of  elemental  carbon  must  be  conserved,  chemistry  ; 
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will  only  effect  the  ratio  of  CO  and  CO^  and  not  the  sum  of  the  two.  As  can 
be  seen  this  ratio  does  not  change  at  the  higher  altitude  and  thus  chemical 
reaction  is  unimportant. 

This  situation  is  reversed  at  the  altitude  of  46  km  where  chemical 
kinetics  does  play  a role.  As  can  be  seen  from  the  histories  in  Fig.  2, 

CO^  is  being  dissociated  to  CO  in  the  near  wake.  This  trend  reverses  at  a 
distance  of  about  80  meters  where  it  can  be  seen  that  the  [CO2I  mole 
fraction  begins  to  increase  while  that  for  [CO]  decreases.  In  general, 
this  byplay  is  a result  of  the  reaction 

CO  + OH  CO^  + H • (7) 

At  the  near  wake  centerline  temperatures  of'-2500°K  the  reverse  of  this 
reaction  is  favored  and  thus  CO^  is  converted  to  CO.  Further  downstream 
in  the  wake,  where  the  temperature  is  somewhat  lower,  the  reaction 
proceeds  in  the  forward  direction  converting  CO  to  CO^.  The  far  wake 
behaviour,  which  is  not  shown,  is  that  the  CO  is  completely  oxidized  to 

CO3. 

This  effect  is  best  seen  by  considering  the  radial  profiles  of  CO  and 
CO^  shown  in  Figs.  3-5  for  wake  distances  of  20,  60  and  100  meters  and 
an  altitude  of  46  km.  An  inert  species,  not  present  in  the  ambient  atmos- 
phere, would  have  a Gaussian- like  radial  distribution  much  like  that  shown 
for  CO  in  Figs.  3-5.  The  wake  temperature  profile  would  also  have  a 
similar  shape.  The  CO^  profiles  shown  in  Figs.  3-5  are  quite  non- 
Gaussian  because  of  chemical  effects.  As  can  be  seen  the  profiles  exhibit 
an  off-axis  peak  which  tends  to  move  toward  wake  center  with  increasing 
distance  behind  the  wake.  This  behaviour  results  from  the  fact  that  the 
forward  direction  of  reaction  (7)  is  favored  in  the  lower  temperature  regions 
of  the  wake. 
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Fig.  5 CO  and  CO2  Mole  Fraction  vs.  Normalized  Wake  Radius,  X = 100 M 
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In  summary,  at  the  highest  altitudes  considered  chemical  effects 
are  unimportant  except  for  the  production  of  trace  species.  On  the  other 
hand  at  the  lower  altitudes  the  chemical  processes  tend  to  convert  free 
radicals  to  more  stable  molecules.  This  conversion  can  be  deterred  in 
laminar  wakes  since  the  temperatures  can  be  large  enough  to  allow  signifi- 
cant equilibrium  concentrations  of  free  radicals.  This,  however,  does  not 
occur  in  the  turbulent  wake  at  the  lowest  altitude  considered  where  the 
temperature  drops  rapidly. 

3.  Vibrational  kinetics 

The  detailed  vibrational  kinetic  scheme  employed  in  the  wake 
calculations  involves  some  117  reactions,  as  seen  in  Appendix  A.  However 
the  general  vibro- relaxation  phenomenology  may  best  be  understood  by 
reference  to  the  vibrational  mode  energy  level  diagram  shown  in  Fig.  h, 
which  was  excerpted  from  Ref.  6. 

As  can  be  seen,  the  CO,  molecule  has  three  vibrational  modes. 
These  modes  are  coupled  to  each  other  by  kinetic  mechanisms  and  the 


mode  is  strongly  coupled  to  the  translational  temperature  because  of 
its  low  vibrational  energy.  CO,  can  vibrationally  communicate  with  other 
molecules,  predominantly  through  the 'j  ^ mode  which  is  in  near  energy 
resonance  with  both  CO  and  N,.  The  v mode  also  has  a very  short 
radiative  lifetime,  — 2.  S msec,  and  this  will  play  an  important  role  in  the 
calculations  as  will  be  seen. 

The  1^2®  molecule  also  has  three  modes  which  are  coupled  together, 
at  low  temperatures  even  more  strongly  than  the  n'.odes  of  CO,.  The 
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Fig.  6 Vibrational  Energy  Spacings  for  Molecules  of  Interest. 
Figure  taken  from  Ref.  6. 


it  can  undergo  a rapid  vibrational  exchange  reaction  with  both  H^O  and 
CO^  and  thus  effectively  provides  a coupling  mechanism  between  these  two 
molecules.  Furthermore,  since  the  radiative  lifetime  of  HF  is  considerably 
shorter  than  that  for  H O,  the  strong  vibrational  coupling  between  H O and 
HF  can  lead  to  vibrational  depopulation  of  H^O  via  HF  radiative  decay. 

Two  simplifying  approximations  have  been  made  in  the  vibrational 
kinetic  scheme  listed  in  Appendix  A.  (1),  In  actuality  the  CO^  >Tiode  is 
degenerate  with  degeneracy  (v  + 1)  in  vibrational  level  v.  However,  the 
vibrational  coupling  reactions  can  be  greatly  simplified  if  it  is  assumed 
that  each  level  of  the  mode  is  doubly  degenerate.  This  assumption  has 
been  made  and  it  can  readily  be  shown  that  the  total  vibrational  energy  in 
is  the  same  in  either  case.  (2),  The  radiative  lifetimes  of  the  H^O 
and  modes  differ  by  a factor  of  twenty  with  that  for  the  mode  being 
the  shorter.  It  is  felt  on  theoretical  grounds  that  the  modes  will  be  strongly 
collis  ionally  coupled  to  each  other.  Therefore  the  lifetimes  of  the  two  modes 
have  been  set  equal  to  each  other  to  ensure  that  decoupling  will  not  occur  in 
the  computer  predictions. 

Shown  in  Figs.  7-11  are  predicted  wake  centerline  temperature 
histories  for  the  altitudes  of  69,  61,  53,  46  and  38  km,  respectively.  In 
each  figure  the  translational  temperature  and  the  vibrational  temperatures 
of  CO,  HF,  ('^2^'  ^^2  ^^3^  displayed.  In  can  be  seen 

that  at  the  highest  altitude  the  vibrational  modes  are  severely  collision 
limited,  i.  e.  , radiative  decay  is  iTiore  rapid  than  collision  deactivation,  so 
that  the  characteristic  vibrational  temperatures  are  significantly  lower  than 
the  translational  temperature.  This  effect  becomes  less  pronounced  with 
decreasing  altitude  and  by  46  km  the  vibrational  temperatures  are  almost 
in  equilibrium  with  translation.  Although  not  shown,  this  same  behaviour, 
ic.  , collision  limiting,  is  also  observed  off  axis  in  the  wake  at  these  altitudes. 
The  vibrational  temperature  can  be  higher  than  translation  in  the  far  wings 
of  the  wake,  however,  the  contaminant  species  concentrations  are  very  low 
there. 
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Fig.  9 Predicted  Wake  Centerline  Temperatures  vs.  Distance,  H = 53  km 
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Fifif.  10  Predicted  Wake  Centerline  Temperatures  vs.  Distance,  H = 46  km 


44 


2b. 00  40.00  60.00  80.00  100.00  p20.00 

fiXIRL  DISTANCE  - (M) 


Fig.  11  Predicted  Wake  Centerline  Temperatures  vs.  Distance,  H = 38km 


Note  the  dramatic  effect  of  the  turbulent  wake. 
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This  behavior  can  be  readily  understood  in  terms  of  the  radiative 
lifetimes  and  kinetics  described  previously.  For  example,  considering 
Fig.  7,  corresponding  to  the  highest  altitude,  the  CO^  (v^)  mode  which  has 
the  shortest  radiative  lifetime  exhibits  the  most  precipitous  drop.  Although 
the  CO2  (v^)  mode  has  a lifetime  of  0.3  sec.,  its  characteristic  temperature 
also  falls  below  the  translational  temperature  since  the  CO^  (Vj)  and  (V2) 
modes  are  strongly  coupled  at  these  wake  temperatures.  HF  also  falls  off 
quite  rapidly  in  line  with  its  short  radiative  lifetime.  The  CO  profile  de- 
creases more  slowly,  although  still  falling  significantly  below  the  translational 
temperature.  This  behavior  is  due  both  to  radiative  decay  and  the  vibrational 
coupling  between  CO  and  CO^  (v^) 

An  effect  of  the  chemical  kinetics  may  be  observed  in  the  H2O  (v^) 
profiles  shown  in  Figs.  7-8.  The  radiative  decay  time  for  the  H^O  (v^) 
state  is  similar  to  that  for  CO  and  furthermore  the  state  can  vibratior.illy 
exchange  with  HF;  thus  one  would  expect  it  to  be  collision  limited.  None- 
theless, the  H2O  (Vj)  temperature  closely  follows  the  translational  tempera- 
ture, It  will  be  recalled  that  H2O  is  a trace  species  at  the  higher  altitudes  and 
that  its  concentration  increases  with  wake  distance  through  chemical  produc- 
tion. Since  the  model  assumes  that  H2O  produced  by  chemical  reaction  is 
created  in  vibrational  equilibrium  with  translation,  the  chemical  reaction 
acts  to  drive  the  H2O  (v^)  vibrational  temperature  towaros  the  translational 
temperature  in  opposition  to  the  radiative/vibrational  kinetics.  Similar 
phenomena,  although  less  pronounced,  affect  the  CO  and  CO2  (V3) 
temperatures  in  the  mid-range  of  altitudes  shown. 

The  temperature  histories  for  the  turbulent  wake  shown  in  Fig.  11 
are  quite  different  from  those  shown  earlier  for  laminar  wakes.  In  this 
case  the  wake  temperature  decays  very  rapidly  because  of  cold  gas  entrain- 
ment. The  vibrational  temperatures  initially  follow  this  rapid  decrease 
but  then  begin  to  freeze  out  at  a wake  temperature  of  approximately  1000°K. 
Although  the  CO2  (V2)  and  H2O  (v^)  temperatures  remain  close  to  the  trans- 
lational temperature  throughout,  the  vibrational  temperature  of  HF,  CO 
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and  (v^)  are  as  much  as  a factor  of  two  higher  than  the  translational 
temperature  at  a distance  100  meters  downstream  of  the  vehicle.  Basically 
what  occurs  is  that  the  ablation  product  species  become  so  diluted  in  the 
wake  that  the  only  important  collision  partner  is  nitrogen.  The  diatomic 
molecules  cannot  be  rapidly  deactivated  by  nitrogen  near  room  temperature 
and  furthermore,  the  coupling  between  the  and  modes  of  CO^,  with 
nitrogen  as  the  collisional  partner,  is  an  inefficient  process  at  these  low 
temperatures. 

The  behaviour  exhibited  in  Figs.  7-11  is  typical  of  polar  molecules 
in  general;  (i.  e.  , molecules  which  have  a dipole  moment  and  so  radiate 
spontaneously,)  and  can  be  readily  understood  from  kinetic  considerations. 
The  characteristic  time  for  a gas  kinetic  collision  between  any  two  mole- 
cules is  defined  by 


= z‘^  = (NVO)'^  (8) 

where  Z is  the  collision  frequency,  N the  total  particle  number  density.  V 

the  mean  random  molecular  speed,  and  a the  gas  kinetic  molecular  cross 

section.  Even  at  the  laminar  wake  temperatures  of  2500°K,  it  is  rare  to 

find  a vibrational  process  more  efficient  than  a tenth  gas  kinetic,  and 

furthermore  such  efficient  reactions  generally  involve  contaminant  species 

having  mole  fractions  of  order  one  percent.  Thus,  althOu^n  collision  ti  ues 

will  vairy  from  reaction  to  reaction  a characteristic  rollisiun  time  for 

vibrational  exchange  or  do  ' i;jion  i'  re-entry  wakes  at  temperatures 

near  2500°K  wo\iU  be  ~ loS, . 

k 

Shown  plotted  in  Fig.  12  are  several  multiples  of  the  gas  kinetic 
collision  time  vs.  altitude  for  a temperature  of  2500°K.  Also  shown  as 
horizontal  lines  are  the  radiative  lifetimes  of  several  of  the  vibratiotuil 
species  of  interest  in  the  present  calculations.  If  the  radiative  lifetime  of 
a species  is  of  the  same  order  of  magnitude  as  the  characteristic  time  for 
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vibrational  exchange/deactivation  collisions,  that  species  will  be  collision 


limited  with  a vibrational  temperature  less  than  that  of  translation.  Thus 


it  can  be  seen  from  Fig.  that  if  the  characteristic  time  for  vibrational 


collision  is  of  order  10  T both  HF  and  the  CO  (v  ) mode  will  be  collision- 

K £i  J 


ally  limited  at  altitudes  above  48  km  as  was  predicted  by  the  detailed 
calculations.  Thus  other  species  which  might  be  introducted  into  the  wake, 
with  similar  radiative  lifetimes,  could  also  be  expected  to  be  collision 


limited. 


The  above  discussion  is  valid  for  the  case  of  a relatively  constant 
temperature  laminar  wake.  In  the  case  of  a turbulent  wake,  with  a rapidly 
decaying  translational  temperature,  the  phenomenon  of  vibrational  freezing, 
rather  than  collisional  limiting,  becomes  of  interest.  Vibrational  freezing 


will  occur  when  there  are  insufficient  vibrational  deactivation  collisions 


within  the  flow  time  to  allow  the  vibrational  temperature  to  drop  with  the 
translational  temperature.  As  a specific  example  consider  the  calculation 


at  38  km  shown  in  Fig.  11.  The  wake  temperature  fell  below  1000  K within 
the  first  twenty  meters  of  wake,  corresponding  to  a flow  time  of  ~ 3 x 10  ^ 


sec.  From  examination  of  Fig.  12  it  can  be  seen  that  this  flow  time 


4 -1/2 

corresponds  to  a collision  time  of  ~ 10  T (note  T scales  as  T ).  Thus, 

In  In  ^ 

*4 

if  the  vibrational  deactivation  efficiency  of  a species  is  £ O (10  ) at  tempera- 


listed in  Appendix  A.  (Note  both  CO  and  can  vibrationally  exchange 


with  N_  however  the  N_  vibrational  temperature  is  also  frozen  at  this  alti- 


tude). 


4.  Radiation  predictions 


Given  the  defined  axial  and  radial  wake  temperature  and  species 


profiles,  infrared  radiation  predictions  may  be  specified  in  any  desired 
bandwidth.  Since  wake  radiation  levels  in  atmospheric  windows  are  of 


tures  below  1000  K,  the  species  will  vibrationally  freeze.  This  was  indeed  the 


case  with  CO,  HF  and  the  CO^  (v^)  mode  as  can  be  seen  from  the  rate  constants 


I 


t 
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particular  interest,  radiation  predictions  have  been  provided  for  the 
respresentative  bandwidths  of  3 - 4|inn,  3-5|-im,  10.  2-12.6^m,  and  17.  2 - 
22.8(im.  The  important  radiators  in  each  of  these  bandpasses  are  listed 
in  Table  II.  This  tabulation  does  not  imply  that  the  nnain  portion  of  the 
bands  listed  fall  in  the  specified  bandwidth:  indeed  this  is  generally  not 
the  case.  For  example,  the  H^O  (v and  CO^  N ^ hands  peak  at  2.  7nm, 
and  the  center  of  the  HF  fundamental  vibration  band  is  at  2.  5^m. 

It  is  to  be  noted  that  the  radiation  processes  in  the  3-  Slim  region 
are  dominantly  vibrational  bands  and  thus  the  predictions  in  this  wavelength 
region  will  be  affected  by  the  vibrational  relaxation  phenomena  discussed 
earlier.  This  is  not  the  case  in  the  higher  wavelength  region  which  is 
dominated  by  radiation  from  pure  rotational  bands.  The  wake  chemical 
kinetics  will  be  most  important  in  this  latter  region  as  well  as  the  fluid 
dynamics  which  establishes  the  translational/ rotational  temperature  of 
the  wake  gases. 

Predictions  of  wake  radiation  per  unit  length  for  these  bandwidths 
at  altitudes  of  69  and  46  km  respectively  are  shown  in  Figs.  13-  20.  The 
history  of  each  band  is  shown  separately.  The  69  km  case  is  characterized 
by  severely  collision  limited  vibrational  temperatures  and  near  frozen 
chemistry  whereas  the  46  km  predictions  correspond  to  near  vibration/ 
translation  equilibrium  and  rapid  chemical  reaction. 

The  69  km  case  corresponded  to  a negligible  H^O  concentration 
relative  to  HF.  Thus  HF  is  the  dominant  radiator  in  the  3-4,  10.  2 - 12.  6 
and  17.  2 - 22.  8iim  bands  as  seen  in  Figs.  13,  15,  16.  The  decrease  in  the  HF 
fundamental  vibrational  band  radiation  with  distance  is  quite  sharp  (see 
Fig.  13)  due  both  to  collis'onal  limiting  and  to  the  fact  that  only  a portion 
of  the  P branch  of  the  band  falls  in  the  3 - 4lim  region.  It  can  be  seen  in 
the  3-  Slim  region,  shown  in  Fig.  14,  that  radiation  from  CO  and 
strongly  dominate  that  from  HF.  This  is  true  at  all  altitudes  and  indeed  the 
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Fig.  13  Wake  Radiation  vs.  Distance  for  the  3-4iim  Bandwidth,  H = 69km 
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Fig.  16  Wake  Radiation  v».  Distance  for  the  l7,2-22.8nm  Bandwidth, 
H = 69  km 
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Fig.  17  SVake  Radiation  vs.  Distance  for  the  Bandwidth,  H =46  km 
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Fig.  19  Wake  Radiation  vs.  Distance  for  the  10.2  - 12.  6 Mm  Bandwidth, 
H = 46  krr. 
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Fig.  20  Wake  Radiation  vs.  Distance  for  the  17.2  - 22.8^m  Bandwidth, 
H = 46  km 
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total  radiation  in  the  3 - 5 vim  bandwidth  is  typically  one  to  two  orders  of 
magnitude  higher  than  that  in  the  3-4vim  bandwidth.  Lastly,  it  should  be 
noticed  that  radiation  from  OH  is  unimportant  in  all  bandwidths. 

The  situation  is  somewhat  different  at  46  km.  It  can  be  seen  in 
Fig.  17  that  the  H^O  (v^)  band  now  dominates  HF  in  the  3 - 4vim  band. 
Furthermore,  the  and  HF  rotational  bands  are  of  about  equal  importance 
in  the  10.  2 - 12.  6 and  17.  2 - 22.  8vim  band  widths.  CO  and  CO^  are  still  the 
dominant  radiators  in  the  3 - Svim  band  and  it  is  to  be  noted  that  the  radia- 
tion histories  in  the  3-4  and  3-  5|lm  bandwidths  decrease  only  slightly 
with  distance  mirroring  the  close  coupling  between  vibration  and  translation. 

The  laminar  wake  radiation  histories  vary  uniformly  between  these 
two  altitudes.  The  turbulent  wake  radiation  histories  which  have  not  been 
shown  exhibit  an  initial  precipitous  decrease,  resulting  from  the  sharp 
near  wake  temperature  drop,  followed  by  a low  level  plateau.  The  net 
effect  of  this  behaviour  is  that  the  total  radiation  in  the  first  100  meters 
of  wake  at  38  km  is  less  than  that  at  46  km.  A more  detailed  discussion 
of  wake  radiation  predictions  for  specific  re-entry  vehicles  of  interest  may 
be  found  in  a companion  report. 
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III.  PARTICULATE  RADIATION 


A.  Introduction 

In  certain  circumstances  it  may  be  advantageous  to  int  oduce  foreign 
substances  into  the  wake  of  a re-entry  vehicle  in  order  to  enhance  the  wake 
infrared  radiation  either  in  a specific  bandpass  or  over  the  whole  of  the 
infrared  region.  Two  possible  techniques  to  achieve  this  goal  would  be 
(1)  the  injection  of  foreign  gases,  and  (2)  the  introduction  of  particulate 
matter  into  the  wake.  Much  of  the  phenomenology  to  be  expected  in  con- 
junction with  foreign  gas  injection  may  be  deduced  from  the  previous  dis- 
cussion on  "baseline"  wake  radiation.  In  particular,  it  appears  that  the 
potential  radiation  available  from  molecules  at  high  altitudes  can  be 
severely  diminished  by  the  effects  of  collision  limiting.  Specifically, 
a criterion  for  optimal  injectant  molecules  would  be  that  they  be  strong 
radiators,  but  conversely,  by  definition, strong  radiators  have  short 
radiative  lifetimes  and  thus  would  be  mcst  affected  by  collision  limiting. 

It  is  the  purpose  of  this  section  to  examine  the  infrared  performance 
of  particles  injected  into  the  wake  at  high  altitudes.  Of  course  it  must 
first  be  ascertained  as  to  whether  one  can  achieve  the  same  radiation 
per  unit  mass  from  particles  as  from  gases.  Be  it  particle  or  gas, 
optimum  radiation  per  unit  mass  is  obtained  under  optically  thin  condi- 
tions. The  absorptivity  of  small  optically  thin  particles  is  generally 
specified  in  terms  of  q,  the  volume  absorption  coefficient,  which  may  be 
related  to  the  absorption  cross  section  by  the  relationship 

* 3 

o = 4/3  ttR  a (9) 

a p 

1 

II 

' where  Rp  is  the  particle  radius.  In  general  this  quantity  a must  be  evaluated 

through  use  of  Mie  theory  and  is  a function  of  particle  size,  wavelength  and 

! 
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4 


the  materials'  index  of  refraction.  The  total  radiation  from  an  optically 
thin  particle  at  a given  frequency  v is 


I = 4/3  TT  R ^ 
P P 


-hv/kT 

a (2  hv  / c ) e 


(10) 


This  is  to  be  compared  to  the  radia  tion  from  a molecule  having  a band 

-1  -2 

centered  at  frequency  v with  bands trength  S in  ama  -cm 


S 

N A V 
o 


2 hv^ 

2 

c 


hv/kT  -1 
(e  -1) 


(11) 


where  is  the  bandpass  of  the  band  and  Nq  is  Loschmidt's  number. 

Both  of  these  quantities  can  be  converted  to  radiation  per  unit  mass, 
J,  upon  introduction  of  the  particles  specific  density  and  the  molecules 
molecular  weight  M,  It  can  then  be  shown  that 


g ^ V M 


2.2x10  p S 
P 


(1 


-hv/kT 

e ) 

1 


(12) 


As  a specific  numerical  example  take  = 250  cm  , M = 100  gms/mole, 
Pp  = 3 gm/cc,  v=  1000  cm' ai 
IR  and  a laminar  wake.  Thence, 


Pp  = 3 gm/cc,  V = 1000  cm'^,  and  T = 2 500°K,  values  typical  of  the  mid 


J /J..  = a / (4S  ) (13) 

p M 

For  the  proper  choice  of  particle  size,  values  of  g of  order  10^  cm"^  may 

be  readily  obtained  and  furthermore  molecular  bandstrengths  of  several 
- 1 -2 

thousand  ama  -cm  are  quite  large.  Thus,  in  principle,  particles  may 
be  used  in  place  of  molecules  as  wake  infrared  sources  with  no  great  increase 
in  weight  penalty. 
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The  remainder  of  this  section  is  concerned  with  the  dynamic  and 
radiative  phenomenology  of  particle/wake  interactions.  The  theoretical 
analysis  of  the  particle/wake  interaction  is  developed  in  section  B and 
applied  to  various  scenarios  in  section  C. 

B.  Flow  Modeling 

The  scenario  considered  is  that  of  small  particles  ejected  from  the 

base  of  a re-entry  vehicle  and  having  relatively  little  or  no  initial  velocity 

relative  to  the  vehicle  (thus  having  a high  velocity  relative  to  ground).  Since 

the  near  wake  velocity  is  low  relative  to  the  vehicle  velocity,  ~0.  1 U , these 

00 

particles  initially  have  significantly  higher  velocity  than  the  wake  molecules. 
Thus,  particle/molecule  collisions  result  in  momentum  and  energy  exchange 
which  have  the  net  effect  of  slowing  and  heating  the  particles. 

The  first  step  in  the  fluid  dynamical  modeling  is  to  define  the  appro- 
priate flow  regime  of  the  particle-gas  interaction.  The  flow  will  be  free 
molecular  if  the  Knudsen  number,  i.e.,  the  ratio  of  the  gaseous  mean  free 
path  to  particle  size,  is  greater  than  ten.  Probstein^^  has  shown  that  in  the 
case  of  hypersonic  flow  the  hyperthermal  rather  than  ambient  mean  free  path 
should  be  applied.  The  hyperthermal  mean  free  path  is  essentially  that 

between  ambient  molecules  and  those  which  have  undergone  perfectly  accom- 

12 

modating  collisions  with  the  particle  and  may  be  defined  by 


X 


H 


V(Tp) 


U 


REL 


X 

a 


(14) 


where  is  the  relative  velocity  between  the  particle  and  the  ambient 

gas,  X^  is  the  ambient  mean  free  path  and  V (Tp)  is  the  mean  Maxwellian 
molecular  velocity  at  the  particle  temperature  T^.  The  hyperthermal  mean 
free  path  is  tabulated  in  Table  III  for  typical  initial  values  of  particle  tempera- 
ture and  relative  velocity.  As  the  particles  slow  down  and  heat  up,  the  hyper- 
thermal mean  free  path  will  approach  that  of  the  wake  which  is  also  tabulated. 
It  is  clear  from  Table  III  that  the  flow  will  be  free  molecular  over  the  full 
altitude  range  of  interest  for  particle  sizes  less  than  1 vini.  Since,  as  will  be 
seen,  the  particle  sizes  required  for  optically  thin  radiation  are  of  the  order 
of  1 yim  or  less,  the  remaining  analysis  will  be  limited  to  the  free  molecular 
flow  regime. 
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TABLE  III 


AMBIENT  AND  HYPERTHERMAL  MEAN  FREE 
PATHS  VS.  ALTITUDE 

T = 2 500°K,  T = 30  0 °K,  = 6 x 10^  c m / s 

p R E L 


H , km 


X,.cm 

a 


X » cm 


69 
6 1 
5 3 
46 


1.1  X 10 

3. 0 X 10'* 

1.2  X 10 

4.4  X 10 


8 . 8 X 10 


2 . 4 X 10 


- 2 


9 . 6 X 10 


3 . 5 X 10 


- 3 


- 3 


3 8 


1.5  X 10 


- 2 


1 . 2 X 10 


- 3 


c c 
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In  the  free  molecular  flow  limit  the  equations  governing  the  heating  and 
deceleration  of  the  particles  may  be  deduced  straightforwardly  from  the 
equations  of  conservation  of  energy  and  momentum  which  describe  the 
"hard  sphere"  collision  between  the  particle  and  a molecule.  The  one 
complication  in  this  analysis  is  the  introduction  of  the  accommodation 
coefficient  a^.  If  a collision  is  perfectly  inelastic,  that  is  if  after  the 
collision  the  gas  molecule  has  the  temperature  and  velocity  of  the 
particle,  the  accommodation  coefficient  is  unity.  At  the  other  extreme, 
a completely  elastic  collision  where  after  the  collision  the  molecule 
maintains  its  original  temperature  and  absolute  velocity, the  accommodation 
coefficient  is  zero.  These  two  limits  obviously  provide  for  quite  different 
collisional  momentum  and  energy  exchanges.  Little  information  on  the  mag- 
nitude of  is  available  although  measured  values  for  air/metallic 

interfaces  are  near  unity^^.  The  situation  is  further  complicated  by  the 
fact  that  the  different  degrees  of  freedom  of  the  molecule  may  have  different 
accommodation  coefficients  and  furthermore  that  the  accommodation  coef- 
ficient may  vary  with  the  angle  of  the  collision. 

In  the  analysis  to  follow  it  will  be  assumed  that  the  accommodation 
efficiency  may  be  defined  by  the  one  coefficient  2^.  The  Schaaf-Chambre 
model  has  been  used  to  describe  the  kinetic  phenomenology.  This  model 
includes  the  effect  of  the  Maxwellian  velocity  distribution  of  the  molecules, 
assumes  the  molecular  component  is  a constant  bath,  and  only  involves 
the  relative  velocity  between  molecules  and  particles.  The  Schaaf-Chambre 
model  has  been  extended  in  the  present  work  to  include  the  two  additional 
effects  of  particle  vapori  ration  and  radiative  loss.  Vaporization  is  also  a 
kinetic  process  in  the  free  molecular  flow  regime  and  the  vaporization 
rate  must  be  defined  in  terms  of  a "sticking"  coefficient  2^,  which  is 
essentially  a measure  of  the  fraction  of  vapor  molecules  which  will  adhere 
to  the  particle  surface  upon  collision. 

If  Q is  defined  as  the  thtermal  heat  content  of  the  particle,  i.e.  , for 
a uniform  temperature  particle 
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Q = 4/3  TT  R ■'  p / C dT 
P P •'  P 


(15) 


where  p and  C are  the  particles  density  and  specific  heal  respectively, 
then  the  particle  heating  phenomenology  is  defined  by  the  relationship 


(16) 


dQ/dt  = dQ/dtj  -dQ/dt.  - dQ/dt  . 

'coll.  I rad  |vap 

where  from  the  Schaaf*Chambre  model 


1 Ml  I 

dQ/dt  ^ A (S)  (2  R T)  4ttR  p„  (T^-T  )a_(vM)/v  (17) 

COM..  p B Ps  ■■  P T 

where  T,  R , p , C and  y are  the  gas  temperature,  specific  gas  constant, 
P P P_ 

density,  specific  hear  at  constant  pressure  and  ratio  of  specific  heats  respect- 
ively, and  S is  the  normalized  relative  velocity  between  the  particle  and 
gas , i.e. , 


S ' (U  - V » / (2  R T) 
P P P 


1 /2 


(18) 


The  quantity  A(S)  is  defined  by 


-1 


A (S)  M8  S)  fS  * Sie  rfc(S)*  1/2  e rf  (S)  1 


(I*?) 


where  erf  is  the  error  function.  is  a modified  recovery  temperature 
defined  by 


T f B (S)  • 1 J T 

r V * 1 


(20) 


and 


B(S)  = 


(2S2  * 1)[S^  ♦Sierfc(S)]  *(2S^-l)erf(S)/2 


(21) 


S**  [S^  -Sierfc  (S)  ♦ erf(S)  / 2 ] 
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(this  last  expression  has  been  corrected  for  typographical  errors  found 

13. 


in  the  original  text  ). 

The  radiative  loss  is  defined  by 


dQ/dt 


RAD 


4tt  R ^ € (n,  R , T ) o T 
P P P P 


(22) 


where  € (n,  R^,  T^)  is  the  size  and  temperature  dependent  particle 


emissivity  which  must  be  deduced  by  use  of  Mie  theory  ( which  requires 
knowledge  of  the  complex  index  of  refraction  n of  the  particle  material) 
and  0 is  the  radiation  constant,  0.57  x 10*  erg/cm'^/sec/  ( K)  . 


Lastly  the  vaporization  loss  is  defined  by 


dQ/dtL,.„  =4ttR  a VAH 

VAP.  P V V V 


(23) 


where 


V = (R  T /2tt) 

V p 


1/2 


(24) 


R is  the  specific  gas  constant  of  the  vapor,  AH  is  the  heat  of  vaporiza 


tion,  and  p is  the  vapor  density  at  the  particle  surface  as  defined  by  the 


Clausius -Clapeyron  relationship,  i.e.  , 


0 = P„/(R  T ) exp  [ah  (T  -T„)/(R  T T_  ) ) 

V B vp  '^‘  vp  B vpB 


(25) 


where  P„  = 10^  dynes /cm^  and  T-  is  the  normal  boiling  temperature  of 
B D 


the  particle  material. 

With  the  definitions  of  Eqs.  (17)  to  (23),  Eq.  (1  1)  simplifies  to 


dT  /dt  = 3(C  0 R )"‘ jA(S)(2R„T'‘^^p  C„  aj  (T  -T  )(y+1)/y  (26) 

P PPP)  R RP^^rp 


g 


e(n,  R,T)aT  -a  Va 

P p p V 


H p } 
V '^vl 


where  it  is  to  be  noted  that,  with  the  exception  of  the  emissivit^',  the  heating 
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T1 


4 


f 

I 
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rate  is  inversely  proportional  to  particle  size. 

The  remaining  equations  required  to  define  the  particle/gas  inter- 
action are  the  momentum  equation 


.1  2 

dV  /dt  = 3 (S)  (8  p R ) p (U  - V ) (27) 

P D P P g g P 

where  is  a modified  drag  coefficient  defined  by  the  Schaaf-Chambre 
model  as 


Cj^(S)  = S‘^  [ (4S"^  + 4S^  - 1)  erf(S)/(2S)  + (2  + 1 ) n 


1/.' 


+ 2/3  S‘^  (ttT 

T p 


(28 


and  the  equation  for  the  rate  of  change  of  particle  size  resulting  from 
vaporization 


d R„/d  t = - a V p /p  . (29) 

p V v p 

It  is  to  be  noted  that  the  phenomenon  modeled  is  only  dependent  upon 
the  relative  velocity  between  particle  and  molecules.  For  simplicity  the 
above  equations  have  been  derived  in  a coordinate  system  where  the  initial 
particle  velocity  is  zero.  This  coordinate  system  is  the  equivalent  to  body 
centered  coordinates  when  the  particles  are  injected  initially  at  rest  rela- 
tive to  the  body.  In  this  case 


U = U - U 

g OD  W 


(30) 


where  U is  the  body  velocity  and  the  wake  velocity,  which  is  typically 

~ 0.1  U in  the  near  laminar  wake.  In  this  instance  time  may  b®  directly 
00 

related  to  distance  behind  the  body  by  the  relationship 
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i 


(31) 


f 


X (t)  = / Vp  dt  . 


If  the  particles  are  injected  with  a finite  velocity  Up  relative  to  the  body 
then  Eqs.  (30)  and  (31)  must  be  replaced  by 


U„  = U - U - U (32) 

g OD  p W 

and  t 

X (t)  = U t + / V dt  . (33) 

P J P 
o 

Eqs.  (26)  • (33)  with  all  peripheral  definitions  have  been  computer 
programmed  and  the  model  has  been  exercised  for  a number  of  scenarios 
as  discussed  in  the  next  section. 
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C,  Particle  Heating  and  Radiation  Predictions 

The  desired  end  result  of  this  analysis  is  a prediction  of  the  near 
wake  infrared  radiation  history  of  specific  particles.  The  choice  of  particle 
material  and  size  will  of  course  be  dictated  by  the  application.  As  alluded 
to  earlier  there  are  two  distinct  forms  of  infrared  wake  enhancement  which 
are  of  interest.  The  first  of  these,  spoofinp,  involves  introducing  an  additive 
which  will  increase  the  IR  radiation  in  one  of  the  atmospheric  windows  but  not 
in  the  others,  for  example,  introducing  a species  which  will  radiate  in  the 
10-12  pm  but  not  in  the  3-5  pm  wavelength  region.  Spoofing  is  a possible 
countermeasure  to  the  use  of  two  channel  infrared  discriminants  as  is  discuised 
in  some  detail  in  the  companion  report^  The  alternate  technique  is  masking 
which  requires  an  additive  which  will  radiate  strongly  throughout  the  entire 
infrared.  The  purpose  of  this  technique  is  to  provide  an  additional  IR  radiation 
source  which  will  exceed  that  of  the  total  uncontaminated  body  and  near  wake 
radiation. 

Particulate  materials  can  readily  be  specified  to  satisfy  either  of 
these  criteria.  For  example,  the  infrared  spectral  volume  absorption  co- 
efficients for  pure  SiO^,  AI2O3  and  carbon  are  shown  in  Fig.  21.  These 
volume  absorption  coefficients  have  been  calculated  in  the  Rayleigh  limit 
corresponding  to  particle  sizes  much  smaller  than  the  radiation  wavelength. 

In  this  limit  a is  related  to  the  materials  complex  index  of  refraction  n by 
the  relationship 


a 


6 TT 


I 


m 


n^  - 1 

2 ^ 
n +2 


(34) 


Detailed  information  on  the  spectral  variation  of  the  complex  index 

of  refraction  of  various  materials  is  generally  hard  to  locate  and  is  clearly 

an  area  which  requires  more  research.  The  values  of  the  complex  index  of 

refraction  used  in  evaluating  Eq.  (34)  were  taken  from  various  sources  and 

do  exhibit  uncertainties  although  the  general  spectral  variations  are  correct. 

14 

The  data  for  AI2O2  was  taken  from  the  review  of  Gal  and  Kirch  and 
are  appropriate  for  a temperature  of  2000°K.  The  data  for  SiO^  is  from 
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the  recent  study  of  Steyer  et  al^^;  the  complex  index  of  refraction  of  this 
material  is  only  available  at  room  temperature.  Lastly,  the  values  for 
carbon  were  taken  from  the  work  of  Stull  and  Platz^^  and  are  appropriate 
for  a temperature  of  2100°K. 

The  main  point  to  be  deduced  from  Fig.  21  is  that  the  materials 
SiO^  and  Al^O^  exhibit  a very  structured  spectra  in  the  infrared,  much  j| 

like  gas  phase  radiators,  whereas  species  such  as  carbon  have  very  smooth 
infrared  spectra.  Thus  materials  such  as  the  former  are  candidates  for  use 
in  spoofing  whereas  materials  such  as  carbon  would  be  used  for  masking 
applications.  (It  should  be  mentioned  that  the  data  shown  in  Fig.  21  are 
for  pure  materials  and  that  impurities  can  act  to  wash  out  much  of  the 
structure  in  the  displayed  spectra).  Since  the  desired  material  properties 
for  masking  are  quite  different  than  those  for  spoofing  the  two  topics  are 
discussed  separately  below. 

1 . Spoofing 

The  literature  was  reviewed  in  order  to  identify  that  material  which 
both  radiates  strongly  in  the  region  of  the  10  - 12.5  ^m  atmospheric  window 
and  very  weakly  in  the  3-5  ^jm  region.  Unfortunately,  the  data  base  is  very 
sparce  and  the  optimal  material  has  not  yet  been  identified.  The  material 
Al^O^  comes  close  to  satisfying  these  requirements  and  has  been  used  as  a 
representative  material  in  the  present  calculations.  The  available  data  on 
the  complex  index  of  refraction  of  Al^O^  has  recently  been  reviewed  by 
Gal  and  Kirch^^  and  is  reasonably  well  defined  for  temperatures  below  the 
melting  point,  2320°K,  although  it  is  quite  uncertain  at  higher  temperatures. 

The  predicted  spectral  variation  of  the  volume  absorption  coefficient 
of  A1_0  in  the  region  of  the  10  - 12.5  unn  atmospheric  window  is  shown  in 

C»  j 

o 

Fig.  22.  The  predictions  are  for  a temperature  of  2000  K and  particle 

radii  varying  from  0.5-5  ijm  and  were  deduced  from  the  Mie  theory 

14 

absorption  efficiency  calculations  of  Gal  and  Kirch  . These  calculations 
were  performed  at  1 jjm  intervals  and  have  been  arbitrarily  connected  by 
straight  lines  thus  causing  the  jagged  appearance  of  the  curves.  Note  that 


for  particle  rad  i larger  than  1,0  ^m  the  absorption  coefficient  decreases 
approximately  inve rsely  with  particle  size.  This  decrease  is  due  to  optical 
thickness  effects  and  thus  true  volume  absorption  is  only  obtained  with 
AI2O2  for  particle  radii  i 1 jjm. 

This  effect  may  also  be  seen  in  the  temperature  dependent  emissivity 

of  Al^O^  shown  for  particle  sizes  between  0.  5 - 5 vim  in  Fig.  23.  These 

1 4 

predictions  were  also  taken  from  the  work  of  Gal  and  Kirch  . The  dashed 
line  intercepts  to  the  predictions  for  radii  of  0.  5 |jm  and  5.  0 nm  represent 
the  extent  of  the  previously  alluded  to  uncertainty  in  the  i'ldex  of  refraction 
at  temperatures  above  2000°K.  It  is  to  be  noted  that  the  emissivities  of  2 
and  Sum  radius  particles  are  approximately  equal,  a behaviour  appropriate 
for  "gray"  bodies.  For  1 uni  and  below  however  the  emissivity  scales 
approximately  inversely  with  particle  radius  as  would  be  expected  for  a 
volume  absorber. 

Although  the  Al^O^  absorption  coefficient  does  not  peak  in  the 
10  - 12,5  uni  region,  see  Fig,  22,  it  is  still  fairly  large  in  the  region  of 
10-11  um  and  thus  can  provide  significant  radiation  in  the  desired  band- 
pass. From  Figs.  22  and  23  it  may  be  deduced  that  optimum  radiation/unit 
mass  may  be  obtained  for  piarticle  radii  1 um  or  smaller.  It  will  be  recalled 
from  Table  III  that  the  flow  of  1 um  particles  will  be  in  the  free  molecular 
regime  at  altitudes  > 38  km, 

A series  of  predictions  for  the  heating  and  radiation  histories  of 

1 um  radius  alumina  particles  have  been  developed  over  the  altitude  range 

of  69  - 38  km.  The  scenario  considered  was  that  of  cold  particles  (T  = 300°K) 

injected  initially  at  rest  relative  to  the  body  into  a constant  temperature  and 

velocity  laminar  wake.  The  gas  properties  were  taken  to  be  that  of 

heated  air  of  temperature  2500°K,  A nominal  reentry  velocity  of  = 6,  6 km/sec 

was  chosen  which,  assuming  that  <«*  0,  1 U^,  corresponds  to  an  initial 

body  centered  gas  velocity  of  U » 6 km /sec.  The  thermodynamic  prop- 

® 17 

erties  of  alumina  were  taken  from  the  JANAF  tables  , The  heat  of  vapori- 
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ization  of  A1  O was  not  readily  available  and  was  estimated  from  Trouton's 

18^^  o 

rule  to  be  860  cal/gm  for  a boiling  temperature  of  3600  K.  One  feature 

peculiar  to  alumina  is  a very  high  enthalpy  of  melting,  ~280  cal/gm.  Since 

the  intention  of  the  calculations  is  to  provide  an  estimate  of  the  generic 

beliaviour  of  typical  materials  the  heat  of  melting  was  taken  to  be  zero  in 

the  present  calculations. 

Predicted  particle  temperature  histories  for  the  altitudes  of  69,  53 
and  38  km  are  displayed  in  Figs.  24  - 26  respectively.  These  predictions 
)iave  been  performed  under  the  assumption  that  the  thermal  accommodation 
coefficient  a>j.  is  unity.  At  all  three  altitudes  three  distinct  cases  were 
considered  in  order  to  delineate  the  importance  of  the  various  losses.  In 
each  figure  the  highest  temperature  curve  corresponds  to  the  case  of  no 
radiation  or  vaporization  losses,  the  next  lowest  curve  to  radiative  loss 
only,  employing  the  Al^O^  emissivity  shown  in  Fig,  22,  and  lastly,  the 
lowest  temperature  curve  corresponds  to  the  inclusion  of  both  vaporization 
and  radiation  losses  with  the  sticking  coefficient  taken  to  be  unity.  The 
particle  velocity  scale  shown  at  the  top  of  each  figure  is  specific  to  the  case 
of  both  radiation  and  vaporization  losses.  However,  the  velocity  histories 
of  the  other  two  cases  are  not  significantly  different. 

Examining  first  Fig.  24,  corresponding  to  an  altitude  of  69  km,  it 

can  be  seen  that  the  initial  particle  heating  is  approximately  the  same  for 

all  three  cases  with  losses  not  becoming  important  until  the  particle 

o 

temperature  reaches  ~ 1800  K,  Furthermore,  if  there  were  no  losses 
the  particle  temperature  would  continue  to  increase  to  temperatures  well 
above  the  wake  temperature.  As  can  be  seen  this  behavioui  is  severely 
affected  by  the  inclusion  of  radiative  losses.  Indeed  the  temperature  profile 
for  this  case  peaks  about  40  meters  downstream  even  though  there  is  still  a 
large  velocity  difference  between  the  gas  and  particle.  The  inclusion  of 
vaporization  decreases  the  temperature  even  more;  so  that  at  this  altitude 
the  peak  particle  temperature  remains  below  the  wake  temperature. 
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In  Fig.  25.  corresponding  to  an  altitude  of  53  km.  the  differences 
between  the  three  curves  are  enhanced.  Because  of  the  increased  gas 
density  the  velocity  difference  between  particle  and  gas  diminishes  more 
rapidly  with  distance.  It  can  be  seen  that  even  the  predicted  temperature  ft> 
the  case  of  no  losses  reaches  a peak  value  and  then  decreases.  This 
occurs  because  the  particle  temperature  is  higher  than  that  of  the  wake 
and  thus  the  particles  can  lose  energy  by  increasing  the  molecule  tempera- 
ture in  perfectly  accommodating  collisions.  If  this  calculation  had  been 
carried  to  larger  distances  the  particle  temperature  would  of  course 
asymptote  to  the  wake  temperature. 

It  is  interesting  to  note  that  the  predicted  temperature  history  for 
the  case  of  radiative  losses  only  is  actually  higher  at  this  altitude  than  at 
69  km.  This  seemingly  contradictory  behavior  occurs  because  the 
particle  velocity  equilibrates  more  rapidly  with  decreasing  altitude.  Thus 
a particle  takes  less  time  to  reach  a given  distance  downstream  as  the 
altitude  decreases  and  therefore  has  less  time  to  radiate  the  excess  particle 
heat  energy  away.  This  effect  is  not  so  apparent  w’hen  vapori r.ation  is 
included  since  a small  change  in  temperature  can  produce  a large  change 
in  vapor  pressure  and  thus  in  the  vaporisation  rate.  Indeed  this  latter 
effect  which  is  controlled  by  the  heat  of  vaporisation  and  the  boiling  tempera 
ture  of  the  material,  see  Eqs.  (23)  - (25).  ultimately  provides  the  limitation 
on  maximum  achievable  AI2O3  Particle  temperatures. 

Similar  behaviour  is  showm  in  Fig.  26  for  an  altitude  of  38  km.  An 
additional  prediction  for  a wake  temperature  of  400°K  has  been  provided  at 
this  altitude  as  a simulation  for  possible  particle  behaviour  in  a turbulent 
wake.  Note  that  in  this  case  the  particle  is  rapidly  heated  to  temperatures 
of  ~ 2 500°K  and  then  slowly  decays  to  the  wake  temperature  within  the 
first  one  hundred  meters  of  wake.  The  velocity  scale  shown  on  this  figure 
is  only  appropriate  for  the  laminar  wake  cases  since  the  "turbulent  wake 
prediction  necessarily  corresponds  to  a higher  gas  density  . 
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These  calculations  may  now  be  applied  to  the  prediction  of  particle 
radiation.  (In  performing  such  predictions  it  is  to  be  remembered  that 
the  particles  decrease  in  size  through  vaporization.  For  the  range  of 
altitudes  shown  the  particles  decrease  at  most  to  70%  of  their  original 
size  in  the  first  100  meters  of  wake.)  The  total  particle  radiation  in  the 
10  - 12. 5 urr  bandwidth  for  optically  thin  particles  is  defined  by  the  relation- 
ship 

12.  5 urn 

Ip=4/3nRp^y  aj(X,  Tp)dX  (35) 

10  um 


where 


• 7 i,  2 ,-5 

j = 2 h c X 


-hc/('  k T) 
e 


(36) 


I 

and  where  the  spectral  variation  of  a is  as  shown  in  Fig.  22.  For  purposes 
of  calculation  it  is  convenient  to  express  the  intensity  in  terms  of  an  average 
value  of  a over  the  bandpass,  i.e., 

I 3 

I Ip  = 4/3  nRp  a j (11.25  um,  ) AX  (37) 

where  bandwidth.  It  <'an  be  shown  that  for  Al^O^  at  temperatures 

of  a few  thousand  degrees  the  average  value  of  a satisfying  Eqs.  (35)  - (37) 
is  5000  cm  (Note  however  that  some  93%  of  this  radiation  will  fall 
between  10-11  ym  .) 

The  predicted  10  - 12.5  jjm  particle  radiation  histories  for  the 
altitudes  considered  are  displayed  in  Fig.  27.  The  particle  radiation  is 
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highest  in  the  near  wake  and  decreases  with  increasing  distance.  This 
decrease  is  not  primarily  due  to  decreasing  particle  temperature.  Indeed 
it  can  be  seen  from  Figs.  24  - 26  that  the  particle  temperature  varies  only 
slightly  in  the  first  100  meters  of  wake.  Rather  it  is  due  to  the  decreasing 
particle  size  since  the  radiation  scales  as  the  cube  of  radius.  This  phen- 
omenon of  particle  shrinkage  has  obvious  advantages  for  applications  where 
it  is  desired  to  enhance  the  infrared  radiation  in  the  near  wake  but  not  in 
the  far  wake. 

One  additional  matter  that  can  be  addressed  is  the  particle  radiation 
resulting  from  the  injection  of  a constant  mass  flow  m^  of  particles  into 
the  wake.  This  cannot  be  deduced  by  a simple  scaling  of  the  results  in 
Fig.  27.  Although  the  particle  mass  flux  is  not  conserved  in  the  wake, 

s 

because  of  vaporization,  the  total  number  flux  of  particles  N is  conserved. 

P 

However,  since  the  particle  velocity  varies  with  distance,  the  local  particle 
density  must  also,  i.e. , 

N = N V A (38) 

P P P 

where  A is  the  cross-sectional  area  over  which  the  particles  are  disseminated. 

Thus  for  a constant  mass  flow  injection  rate  the  number  of  radiating  particles 

per  unit  length  of  wake  will  scale  inversely  with  particle  velocity. 

As  a numerical  example  radiation  predictions  have  been  provided  at 

various  altitudes  for  initial  injection  rates  which  produce  a total  irradiation 

of  20  w/sr  in  the  first  30  meters  of  wake.  The  required  m to  satisfy  this 

P 

criterion  is  deduced  from  the  relationship 

/30  m 

I /V  dx  = 20  w/sr  . (39) 

P P 
o 

The  rerulting  radiation  histories  along  with  the  required  injection  mass  flows 
for  1 pm  alumina  particles  are  displayed  in  Fig.  28.  It  can  be  seen  that  the 
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required  mp  scales  inversely  with  altitude  reflecting  the  increased  rates  of 
both  velocity  equilibration  and  vaporization  which  occur  with  decreasing  altitude. 

2.  Masking 

There  has  to  date  been  no  detailed  study  performed  to  determine  which 
particle  properties  such  as  size,  emissivity  and  specific  density,  will  lead  to 
optimum  radiation  per  unit  mass  for  masking  applications.  In  the  present 
analysis  two  specific  particle  types  have  been  arbitrarily  chosen  in  order  to 
examine  the  altitude  dependent  infrared  radiation  histories  of  strongly 
radiating  particles  . The  two  cases  considered  are  (1)  Sum  radius  carbon 
particles  which  exhibit  an  emissivity  of  approximately  unity  throughout  the 
IR  and  (2)  0.55  ^m  radius  zirconia  particles  having  a characteristic  constant 
emissivity  of  0.22  in  the  IR. 

The  thermodynamic  properties  of  these  two  materials  were  taken  from 
17 

the  JANAF  tables  . Carbon  has  a boiling  temperature  (more  properly  sub- 
limation) of  4500°K  and  a heat  of  vaporization  of  171  kcal/mole  uhile  the 
respective  values  for  ZrO^  are  4550°K  and  149  kcal/mole.  It  was  found 
in  the  course  of  the  calculations  that  these  boiling  temperatures  are  suffi- 
ciently high  so  that  vaporization  plays  essentially  no  role  in  the  particle 
heating  phenomenology. 

The  approach  taken  in  performing  the  masking  predictions  was  quite 
distinct  from  that  of  the  spoofing  analysis  of  the  previous  section.  Masking 
is  best  achieved  by  employing  a devico  for  the  generation  of  hot  particles 
so  that  a radiation  signature  can  be  provided  at  very  high  altitudes  (exo- 
atmospheric)  where  heating  by  the  particle/gas  interaction  is  negligible. 

In  this  high  altitude  vacuum  limit  the  total  radiation  attainable  from  the 
particle  is  essentially  equal  to  its  initial  heat  content  at  the  time  of  injection. 
The  present  calculations  were  directed  towards  determining  the  degree  of 
enhancement  of  this  vacuum  limit  radiation  which  can  result  from  the  gas/ 
particle  interaction  under  endo -atmospheric  conditions . To  this  end  calcula- 
tions were  performed  over  the  altitude  range  of  91  - 61  km  under  the  assump- 
tions of  an  initial  particle  temperature  of  1840°K  and  an  initial  particle 
separation  velocity  (from  the  body)  of  0.6  km/sec.  Wake  conditions  were 
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taken  to  be  those  of  2500°K  air  with  = 6 km/sec  (see  Eqs.  (32)  - (33)  ), 

The  predicted  temperature  histories  for  5 vjm  radius  carbon 
particles  at  altitudes  of  91 « 76  and  61  km,  and  in  the  vacuum  limit  are 
shown  in  Fig.  29.  As  can  be  seen  there  is  little  difference  between  the 
predictions  at  91  km  and  the  vacuum  limit  implying  that  the  gas /particle 
interaction  is  unimportant  at  this  altitude.  The  particle  temperatures 
are  enhanced  at  the  lower  altitudes  shown  although  the  particle  temperature 
is  everywhere  lower  than  the  wake  temperature. 

IR  radiation  predictions  for  these  altitudes  for  wavelengths  of  1 , 5 
and  15  um  are  shown  in  Figs.  30  - 32  respectively.  These  predictions  are 
for  the  total  integrated  radiation  at  a given  wavelength  per  unit  particle 


mass  flow  K.  vs.  distance,  where 


P P o 


/V  dx 
X P 


where  € i*  the  emissivity  and  the  black  body  function  (the  units  of 
are  joules /(gm-s r-nm)  ).  As  can  be  seen  the  radiation  enhancement  factor 
is  a strong  function  of  wavelength  being  as  large  as  a factor  of  6 at  1 um 
and  close  to  unity  at  15  ^im.  This  wavelength  dependence  is  primarily  a 
result  of  the  s{>ectral  variation  in  the  black  body  function  at  the  temperatures 
of  interest. 

Predicted  temperature  histories  for  0.  55  um  ZrO^  particles  are 
shown  in  Fig.  33  for  the  same  altitude  range.  The  displayed  zirconia 
temperatures  are  somewhat  higher  than  those  shown  earlier  for  carbon. 
Indeed,  close  to  the  body,  the  particle  temperature  exceeds  the  wake 
temperature  at  61  km  altitude.  The  corresponding  mass  flow  normalized 
radiation  predictions  for  wavelengths  of  1 and  5 um  are  shown  in  Figs.  34 
and  35  respectively.  It  can  be  seen  that  the  radiation  enhancement  factor  is 
considerable  in  this  case,  being  as  much  as  a factor  of  60  at  1 uni. 

Although  not  shown  the  radiation  predictions  at  15  um  exhibit  negligible 
radiation  enhancement  much  like  the  predictions  for  carbon  particles. 
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Fig.  29  Predicted  Carbon  Particle  Temperature  v«.  Distance  at 
Several  Altitudes. 
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Fig.  30  Integrated  Carbon  Particle  Radiation  Per  Unit  Mass  Flow  vs. 
Distance,  X = 1 
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Fig.  33  Predicted  Zirconia  Particle  Temperature  vs.  Distance  at 
Several  Altitudes 
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It  ihould  be  noted  that  these  radiation  predictions  can  be  pro. 
vided  for  any  wavelength  and/or  bandpass.  In  particular  predictions  for 
the  bandpasses  corresponding  to  atmospheric  windows  may  be  of  interest 
for  systems  applications.  As  a specific  example  the  predicted  total 
normalized  ZrO^  particle-  radiation  for  the  3-5  um  wavelength  region  is 
shown  in  Fig.  36  for  both  an  altitude  of  61  km  and  the  vacuum  limit.  The 
radiation  enhancement  factor  fur  this  case  is  approximately  a factor  of  five, 
somewhat  higher  than  that  shown  for  5 um  in  Fig.  35  because  of  the  in- 
creasing radiation  enhancement  with  decreasing  wavelength. 

D.  Discussion  and  Conclusions 

In  summary  it  has  been  demonstrated  in  the  spoofing  analysis  that, 
given  the  proper  choice  of  particle  material  and  size,  significant  infrared 
signatures  can  be  produced  in  the  near  wake  of  a re-entry  vehicle  through 
the  mechanism  of  introducing  cold  particles  into  the  flow.  These  particles 
which  are  initially  injected  at  or  near  the  body  velocity  (in  ground  based 
coordinates)  are  heated  through  decelerating  collisions  with  the  wake 
molecules.  Many  materials  are  available  which  exhibit  very  structured 
molecular -like  spectra  in  the  infrared.  Unfortunately,  the  information 
required  to  define  optimum  materials  for  specific  applications  is  not 
generally  available. 

In  the  case  of  the  masking  predictions  it  has  been  demonstrated 
that  there  can  be  a significant  enhancement  in  the  IR  radiation  signatures 
from  PCM  generated  particles  when  endo-atmospheric  effects  become 
important.  The  level  of  radiation  enhancement  is  both  altitude  and 
wavelength  dependent  and  is  most  prominent  in  the  near  infrared.  The 
particle  properties  required  for  optimal  radiation  enhancement  have  not 
yet  been  defined. 

In  extrapolating  the  results  of  this  modeling  effort  to  particular 
systems  there  is  one  limit,  tion  which  must  be  considered,  dealing  with 
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Fig.  36  Integrated  Zirconia  Particle  Radiation  Per  Unit  Mass  Flow 
vs.  Distance,  X = 3-5iim 
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the  total  particulate  mass  flow  required  to  produce  a desired  radiation 
intensity.  The  model  used  in  this  work  assumes  that  the  wake  may  be 
treated  as  a constant  bath.  However,  the  particles  initially  have  a high 
kinetic  energy  and  it  can  be  shown  that  > 75%  of  this  energy  is  transferred 
to  the  wake  molecules  with  the  remainder  being  dissipated  among  radiation 
and  vaporization  losses  and  the  sensible  heat  content  of  the  particles. 

Thus,  if  the  particle  mass  flow  rate  is  such  that  the  initial  particulate 
enthalpy  flow  is  of  the  same  order  or  greater  than  the  wake  molecular  en- 
thalpy flow  the  particle /gas  interaction  could  produce  a significant  pertur- 
bation on  wake  temperature  and  velocity.  Although  this  effect  would  not 
invalidate  the  basic  heating/radiation  phenomenology  described  in  this 
work  it  could  provide  for  a variation  in  radiation  per  unit  mass  as  a function 
of  rop. 
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APPENDIX  A 

CHEMICAL/VIBRATIONAL  KINETIC  PACKAGE 

The  chemical/vibrational  kinetic  scheme  employed  in  the  calcula- 
tions  is  listed  on  the  following  pages.  With  the  exception  of  OH,  NO  and 
each  of  the  important  vibrational  modes  of  a molecule  were  carried  as 
distinct  species.  The  code  use'*  i i designating  those  species  is  listed  in 
Table  1. 

It  will  be  noted  that  in  several  cases  a specific  reaction  will  be  listed 
twice.  In  these  instances  the  total  rate  constant  fur  the  reaction  is  the  sum 
of  the  two  listed  rate  constants.  All  rate  constants  are  specified  in  units  of 
cm^/sec  or  cm^/sec*^.  The  rates  for  spontaneous  radiation  are  listed  in 
units  of  sec  . 
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TABLE  I 

SPECIES  DESIGNATION  CODE 


NAME 

SPECIES 

1 

CO 

CO  (v  = 0) 

2 

CO* 

CO  (v  = 1 ) 

3 

»2 

o 

II 

<M 

X 

4 

H 

H,  (V  = 1 > 

2 

5 

HF 

HF  (V  = 0) 

6 

HF 

HF  (V  - 1 ) 

7 

^2 

(V  = 0) 

8 

N,  (V  - n 

9 

°2 

O^  (V  : 0) 

10 

°2 

Oz  1) 

1 1 

cooo 

CO^  (000) 

12 

ClOO 

CO^  (100) 

13 

COl  0 

CO^  (010) 

14 

C020 

CO^  (020) 

15 

C030 

CO^  (030) 

16 

cool 

CO^  (001) 

1 7 

HOOO 

H^O  (000) 

18 

HlOO 

H^O  (100) 

19 

HOlO 

H^O  (010) 

20 

H020 

H^O  (020) 

21 

HOOl 

H^O  (001) 

) 
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TABLE  II 


Vibrational  Kinetics  k = 

AX® 

AT  ex 

p (C/T 

+ D/T 

E/T^'^ 

1 

REACTION 

A 

B 

c 

D 

E ' 

1. 

CO  ^ CO  i!  CO  * CO 

6.67 

X 

10-* 

0 

0 

-208.  3 

0 

2. 

Nz  ^ ^ . N2 

1. 00 

X 

10-"> 

0 

0 

0 

0 1 

3. 

CO  * O^  ^ CO  + O2 

6.67 

X 

.0-“ 

0 

0 

-208. 3 

0 

4. 

CO*  « ^ CO  * N’2 

6.67 

X 

10-« 

0 

0 

-208.  3 

° t 

5. 

CO  * H2  ? CO  + 

9.  00 

X 

10-'’ 

0 

0 

-1  18.  8 

0 

6. 

CO*  + HOOO  ;*  CO  * HOOO 

2.  55 

X 

10-“> 

0 

0 

-32.4 

0 

7. 

CO  * H ^ CO  ♦ H 

1.80 

X 

10'^ 

0 

0 

.208.  3 

0 

8. 

CO  * 0 ^ CO  ♦ 0 

2.70 

X 

.0-’ 

0 

0 

.69.  7 

0 , 

9. 

Co'  * F ? CO  * F 

5.05 

X 

io-’ 

0 

0 

-95. 91 

0 1 

10. 

‘ -N, 

8.53 

X 

■ 0-’ 

0 

0 

-273.  1 

0 1 

J 

11. 

Nj  • O3  « Nj  . O2 

8.  53 

X 

.o-’ 

0 

0 

.273.  1 

1 

0 

12. 

• CO  N • CO 

8.  53 

X 

10'^ 

0 

0 

-273.  1 

1 

0 

13. 

• H,  # ■ Hj 

3.  92 

X 

10*^ 

0 

0 

- 164.  2 

0 

-10 

14. 

N • 0 < Nj  • 0 

1.  10 

X 

10 

0 

0 

-69.  Q 

0 

.11 

1 5. 

■'■'2  • F ? Nj  • F 

5.  50 

X 

10 

0 

0 

.69.  9 

0 

16. 

Nj  • HOOO  .•  N,  • HOOO 

1.  10 

X 

I0*‘® 

0 

0 

.69.  9 

0 

17. 

°2‘°2^°2‘°2 

1.20 

X 

10-' 

0 

0 

-133.  1 

( 

0 i 

1 

18. 

0^  . « O2  . Nj 

4.  81 

X 

10*^ 

0 

0 

-169.6 

0 ! 

19. 

0 • CO  S O2  ♦ CO 

4.81 

X 

10-^ 

0 

0 

- 169. 6 

0 i 

20. 

0/  * H2  ;•  0^  . 

1.  30 

X 

10"' 

0 

0 

-76.  8 

0 

21. 

0^  ♦ HOOO  ? O2  * HOOO 

4.  10 

X 

10-'^ 

0 

0 

-76.  8 

0 

22. 

°2  " 0 * 0 

4.  50 

X 

10-15 

0.  5 

0 

0 

0 * 

23. 

0^  * F ^ O2  + F 

2.20 

X 

10-‘5 

0.  5 

0 

0 

° 1 

24. 

«2  * »2  ^ "2  ^ »2 

2.20 

X 

10"' 

0 

0 

-1 18. 9 
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TABLE  II  (Cont'd) 


REACTION 

A 

- 

B 

D 

25. 

+ H ? H2  + H 

1 . 70  X 

lO'*** 

0.  5 

0 

0 

26. 

H/  4 Nz  ^ N2 

8. 60  X 

10-9 

0 

0 

- 148.  4 

27. 

H.,"  * HOOO  ^ H,  * HOOO 

1 . 70  X 

10-‘5 

0.  5 

0 

0 

28. 

HO  1 0 t HOOO  HOOO  4 HOOO 

1 . 60  X 

io-'2 

0.  5 

0 

0 

29. 

H020  4 HOOO  ii  HOlO  * HOOO 

3.20  X 

1 2 

0.  5 

0 

0 

30. 

HOlO  * 0 ;!  HOOO  • 0 

1 . 00  X 

10-" 

0.  5 

0 

0 

31. 

HOlO  ‘ F ? HOOO  ♦ F 

5. 00  X 

10-'^ 

0.  5 

0 

0 

32. 

Hf)20  * 0 HOlO  • 0 

2 . 00  X 

10-" 

0.  5 

0 

0 

33. 

H020  * F HOlO  F 

1 . 00  X 

I0-" 

0.  5 

0 

0 

34. 

HOlO  * H HOOO  * H 

1.00  X 

10-" 

0.  5 

0 

0 

35. 

H020  ^ H ? HOlO  ♦ H 

2 . 00  X 

.0-" 

0.  5 

0 

0 

36. 

HOlO  * N'z  HOOO  + N’2 

4.  80  X 

10-" 

0 

0 

-43.  0 

3 7. 

H020  * ^ HOlO  4 

9.  60  X 

10*'  ' 

0 

0 

-43.  0 

38. 

COlO  4 cooo^cooo*  COOO 

4 . 64  X 

10-'° 

0 

0 

-76.  75 

39. 

C020  4 COOO  COlO  4 COOO  9.  28  X 

10-10 

0 

0 

-76.  75 

40. 

C030  4 COOO^ C020  4 COOO 

1 . 39  X 

10"*^ 

0 

0 

-76.  75 

41. 

COl  0 4 ? COOO  ♦ N2 

6. 69  X 

lO*"' 

0 

0 

-84. 07 

42. 

C020  4 ;♦  COlO  4 

1 . 34  X 

10-° 

0 

0 

-84.  07 

43. 

C030  4 N,  ^ C020  4 N, 

2.01  X 

10-° 

0 

0 

-84.  07 

44. 

COlO  + 0^  ^ COOO  4 

6. 69  X 

io-‘° 

0 

0 

-84. 07 

-9 

45. 

C020  4 ;?  COlO  4 

1 . 34  X 

10 

0 

0 

-84. 07 

46. 

C030  + i?  C020  4 

2.01  X 

10-° 

0 

0 

-84.  07 

47. 

COl  0 4 H000i?C000  4 HOOO 

3.22x 

10-'^ 

0 

0 

22.  91 

48. 

C020  4 H000;»C010  4 HOOO 

6.  44  X 

10-‘  ^ 

0 

0 

22.  Q1 

49. 

C030  4 HOOO?  C020  4 HOOO 

9.  66  X 

10-'  ^ 

0 

0 

22.  Q1 
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TABLE  U (Cont'd) 


REACTION 

A 

— 

B 

c 

D 

E 

SO. 

COlO  * H, COOO  • 

7.  60 

X 

10-*’ 

10-" 

io-“ 

0 

0 

-7.  1 3 

0 

SI. 

C020  ♦ M,  ^ COlO  • 

1.  S2 

X 

0 

0 

.7.  1 3 

0 

S2. 

C030  ♦ Ml  ? C020  • M, 

2.28 

X 

0 

0 

-7.  1 3 

0 

S3. 

COlO  ♦ O ? COOO  • O 

2.  32 

X 

10-'^ 

0 

0 

.76. 7S 

0 

S4. 

COlO  • F ^ COOO  • F 

1 . Itj 

X 

.Q 

10 

0 

0 

-76. 7- 

0 

ss. 

C020  ♦ O ? COlO  • O 

4.  b4 

X 

.Q 

10 

0 

0 

- 7f>.  75 

0 

S6. 

C020  • F ? COlO  • F 

2.  32 

X 

.Q 

10 

0 

0 

-76.  7=< 

0 

S7. 

C030  ♦ O ? C020  • O 

0.  «»6 

X 

lO'*^ 

0 

0 

-76. 

0 

S8. 

C030  • F ? C020  • F 

3.  48 

X 

10*‘* 

0 

0 

-76.  7=^ 

0 

S'!. 

*^010  . CO  ? COOO  • CO 

6.  o'i 

X 

lO-'O 

0 

0 

.84.  07 

0 

<)0. 

C020  ♦ CO  COlO  • CO 

1 . 34 

X 

10 

0 

0 

.84. 07 

0 

61. 

C030  ♦ CO  ? C020  • CO 

2.  01 

X 

.Q 

10 

0 

0 

.84. 07 

0 

62. 

MF  • HF  ? HF  • MF 

1.  30 

X 

10**’ 

0 

0 

0 

0 

6 3. 

HF  • MF  ? MF  • MF 

2. 00 

X 

io-« 

0 

0 

.12  1.81 

0 

64. 

MF  • J?  HF  • 

1.24 

X 

io-“ 

0 

0 

-5  7.  56 

0 

6S. 

MF  • O,  ? HF  ♦ O, 

2.08 

X 

10’*  * 

0 

0 

-68. 02 

0 

66. 

HF  * NO  t HF  • NO 

1.20 

X 

io-'° 

0 

0 

.60.  20 

0 

6 7. 

HF  ‘ F ? MF  • F 

1. 40 

X 

io-« 

0 

0 

-70. 88 

0 

68. 

MOlO  • O,  ? O,  ♦ MOOO 

l.OO 

X 

10-*' 

0 

0 

0 

0 

69. 

cool  * CO  j»  CO  * COOO 

1 . S6 

X 

io-‘* 

.10 

0 

0 

-30.  12 

0 

70. 

N’2  * 02  «?  02  * N'2 

* 

1. 74 

X 

10 

-6 

0 

0 

.124. 00 

0 

71. 

cool  ♦ N j*  N ♦ COOO 

1.71 

X 

10 

-14 

0 

0 

-17S. 30 

0 

72. 

cool  * N^  iJ  N^  * COOO 

6.  07 

X 

10 

0 

0 

1 S.  27 

0 

73. 

N^  * CO  t CO^  ♦ 

1 . 78 

X 

10"^ 

0 

0 

-209.  90 

0 

74. 

N^*  > CO  t CO"^  * 

6.  98 

X 

.0-” 

0 

0 

-25.  60 

0 
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TABLE  11  (Cont'd) 


REACTION 

B 

C 

D 

E 

75. 

CO  ♦ O,  ? O,  • CO 

1.00 

X 

io-‘“ 

0 

0 

.104.  00 

0 

7o. 

• HOOO  ? HOOl  • H, 

1.67 

X 

2.  00 

0 

0 

0 

77. 

H,  • HOOO  iJ  HlOO  • H, 

1.67 

X 

lO-**^ 

2.  00 

0 

0 

n 

78. 

w Z 

O . COOO  O,  • C020 

1.40 

X 

10-*' 

1.  50 

0 

0 

0 

79. 

H^  • MF'  • HF 

1.  50 

X 

10-*^ 

10-'- 

0 

0 

0 

0 

80. 

HF  ‘ COOO  cool  • HF 

1.20 

X 

0 

0 

0 

0 

81. 

HF  - HOOO  HF  • HOOl 

5.  00 

X 

10-" 

0 

0 

0 

0 

82. 

cool  • HOOO  i*  C030  * HOOO  4.  71 

X 

io--*« 

4.  54 

2082. 

4 54.  06  . 

1729.2 

83. 

cool  • H^  C030  • H^ 

4.27 

X 

10*' 

.6.  34 

.3013. 

. 378.  66 

1 4 1 5.  K 

84. 

cool  • j*  C0  30  ♦ 

4.  12 

X 

10^^ 

-5.  89 

-4223. 

.672. o 7 

2682. 9 

85. 

cool  ♦ ? C030  • ©2 

4.  12 

X 

10^^ 

-5.  89 

.4223. 

.672.  67 

2oS2. 9 

86. 

cool  * CO  ? C030  * CO 

6.87 

X 

10^' 

-5.  89 

.4223. 

.672 . 6 7 

2ob2.  ' 

87. 

cool  • COOO  j*C030  ‘ COOO 

'L  16 

X 

10^' 

.5.  89 

.422  3. 

.672. 67 

2<,82. 9 

88. 

cool  ♦ O ? C030  * O 

1.20 

X 

lo-*-* 

0.  50 

0 

0 

n 

89. 

cool  * F j*  C030  • F 

1.90 

X 

10-*' 

0.  50 

0 

0 

0 

90. 

cool  * H ;♦  C030  * H 

6.  00 

X 

10-*^ 

0.  50 

0 

0 

0 

91. 

cool  ^ HF  ? C030  * HF 

4.00 

X 

10-*' 

0 

0 

0 

0 

92. 

HlOO  * HOOOi*  H020  + HOOO 

1.60 

X 

10-*^ 

0.  50 

0 

0 

0 

93. 

HOO  1 ♦ HOOO  ?H020  * HOOO 

l.bO 

X 

10-*2 

0.  50 

0 

0 

0 

94. 

HlOO  * O ^ H020  + O 

1.00 

X 

10-** 

0.  50 

0 

2 

0 

95. 

HlOO  * F H020  » F 

5.00 

X 

10-*^ 

0.  50 

0 

0 

0 

96. 

HOOl  + O ^ H020  + O 

1.00 

X 

1Q-*  * 

0.  50 

0 

0 

0 

97. 

HOOl  * F ^ H020  + F 

5.00 

X 

10-*^ 

0.  50 

0 

0 

0 

98. 

HlOO  + H H020  + H 

2.00 

X 

10"*  * 

0.  50 

0 

0 

0 

99. 

HOOl  + H J*  H020  + H 

2.00 

X 

1Q-*  * 

0.  50 

0 

0 

0 
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TABLE  11  ICont’fB  i 


R EACTION 

A 

B 

c 

1) 

E 

I 00. 

H020 

• MOOOtfMOIO  • MOlO 

1 . 00  X 

,o-‘“ 

O 

• 

o 

0 

0 

0 

101. 

moo 

• ? H020  • N, 

4.  80  X 

10-" 

0 

0 

.4  3.  00 

0 

102. 

moo 

• M ? MOOl  • M 

1. 00  X 

10-‘“ 

0.  SO 

0 

0 

0 

1 OJ. 

HOOl 

• a H020  • 

4.  80  X 

10‘  “ 

0 

0 

o 

o 

0 

I 04. 

C I 00 

• .M  ? C020  • M 

1. 80  X 

10' 

0.  so 

0 

0 

0 

1 OS. 

CO  30 

• C000JC020  • COlO 

00  X 

,o-'> 

0.  so 

0 

0 

0 

1 Oo. 

C020 

• COOOjCOlO  • COlO 

4.  00  X 

,o-'> 

0.  so 

0 

n 

(' 

107. 

C030 

• COl  0JC020  • C020 

1. 20  X 

10'*“ 

0.  '0 

0 

0 

0 

TABLE  m 


Spontancou*  Radiative  Decay 


REIACTION 

1 . Co“  CO  • h>. 

2.  cool  ■’  COOO  • fiv 

3.  COlO  COOO  • h'. 

A.  C020  COlO  • h, 

C030  -•  C020  • hv 
o.  HO  10  HOOO  • h'. 

7.  H020  HOlO  • h- 

S.  moo  HOOO  • h>. 

HOOI  -•  HOOO  • h . 
10.  Hf'  - HF  • h. 


. 1 

A sec 

^ ^ » 

3.  03  X lo“ 

3.<>.’  X 10“ 

1 . S'  X 10 

3.  10  X 1 

4. o ' X 1 

2.24  X 1 

4 . 4 » X i r ’ 

3.  X lo‘ 

3.  -Hi  X 10  ‘ 

> 

2. 00  X 10“ 
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TABLE  IV 


D 


Chemical  Kinetio  k - AT  exp  (C/T) 


REDACTION 

A 

B 

c 

1 . 

O • 

O • O,  ? O,  • O, 

3.0  X 10*^° 

.1.0 

0 

O • 

O - O,  ;•  O,  • O, 

3.00 

X 

10*'° 

-1.0 

.224'' 

3. 

C • 

O • ^2  " °2  ’ ^2 

7.  50 

X 

io-^‘ 

-1.0 

0 

4 , 

O • 

O • N , ^ O,  . N , 

7.  50 

X 

10-'' 

-1.0 

.2240 

S , 

CO 

• OH  i COOO  • M 

^.00 

X 

io-‘“ 

-0.  5 

-300 

CO 

• OH  i cool  • H 

6.  00 

X 

10-'- 

-0.  5 

.3680 

« 

CO 

• OH  ? C 100  • H 

6.00 

X 

10*'’ 

-0.  5 

.2  300 

CO 

• OH  f COlO  • H 

1.20 

X 

10*>‘ 

-0.  5 

.1260 

CO 

• OH  f C020  • H 

1.20 

X 

10-" 

.0.  5 

.2  180 

10. 

CO 

• OH  ? C030  • H 

1.20 

X 

10  - ' ' 

.0.  5 

.3080 

1 1. 

CO 

• OH  ;!  COOO  • H 

6.00 

X 

10*'" 

-0.  5 

.300 

y 

CO 

• OH  jf  COOl  • H 

b.  00 

X 

10*'“ 

-0.  5 

.3680 

. 3. 

CO 

• OH  ? ClOO  • H 

6.  00 

X 

10-'' 

-0.  ^ 

-2  300 

! -4. 

CO 

• OH  r COlO  • H 

1.20 

X 

10-" 

-0.  8 

-12b0 

1 

CO 

• OH  t C020  • H 

1.20 

X 

10  - ' ' 

-0.  5 

.2  180 

\f>. 

CO 

• OH  f C030  • H 

1.20 

X 

lo'" 

-0.  5 

.3080 

1 7. 

o • 

CO  • M ? COOO  • M 

3.00 

X 

10'" 

0 

-2170 

1 8. 

o • 

CO  • M r COOl  • M 

3.00 

X 

10-" 

0 

-2  170 

I 

M • 

H • M ? H,  • M 

2.50 

X 

10-'° 

-1.0 

0 

20. 

H • 

H • V ? H,  • M 

2.50 

X 

10-'° 

-1.  0 

-6000 

2 1. 

M • 

OH  K,  • O 

1.40 

X 

10-'^ 

l.O 

-3500 

2 2 

M • 

OH  r IL,*  • O 

1.40 

X 

10-'-* 

1.0 

.9500 
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TABLE  IV  (Cont'd) 


REACTION 

A 

B 

c 

23. 

+ OH  HOOO  + H 

6.  00  X 

10'** 

0 

-2600 

24. 

H2  + OH  HOlO  + H 

6.  00  X 

10‘*  * 

0 

-4900 

25. 

H^  + OH  H020  + H 

6.  00  X 

10'** 

0 

-9400 

26. 

H^’"  + OH  j!  HOOO  + H 

6.  00  X 

10‘*  * 

u 

-2600 

2 7. 

H^"^  + OH  ;*  HOlO  + H 

6.  00  X 

10'*  * 

0 

-4900 

28. 

H^  ''  + OH  ;♦  H020  + H 

6.  00  X 

10’ 

0 

-9400 

29. 

OH  f OH  ^ HOOO  + 0 

8.  20  X 

'0’*‘ 

.0.  33 

-550 

30. 

OH  + OH  ;!  HOlO  + 0 

8.20  X 

10'*  * 

.0.  33 

-2850 

31. 

OH  + OH  ;*  H020  + 0 

8.20  X 

10‘‘‘ 

.0.  33 

-5000 

32. 

OH  + OH  ;*  HlOO  + 0 

8.20  X 

10*‘‘ 

.0.  33 

.5800 

3 3. 

OH  + OH  ;!  HOOl  + 0 

8.20  X 

10‘*  * 

.0.  33 

-5950 

34. 

H + OH  + M HOOO  + M 

6.  70  X 

10-^^ 

.2.  0 

0 

35. 

H f 0 0 + OH 

3.  70  X 

10-‘° 

0 

-8450 

36. 

H + O^*  0 + OH 

3.  70  X 

lo-*'^ 

0 

-8450 

37. 

F + H^  ^ HF*  + H 

2.  60  X 

10-*^ 

0 

-800 

38. 

H + F^  ^ HF*  + F 

2.20  X 

lO-*'^ 

0 

-1200 

39. 

F^+MJI  F+F  + M 

8.  30  X 

10-** 

0 

-17600 

40. 

HF  + MJI  H + F + M 

1 . 80  X 

10’^ 

.1.0 

-67500 

-88- 


